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We present a detailed analysis of the phase diagram of antiferromagnets with competing exchange-driven 
and field-induced order parameters. By using the quasi- ID antiferromagnet BaCu2Si207 as a test case, we 
demonstrate that a model based on a Ginzburg-Landau type of approach provides an adequate description of 
both the magnetization process and of the phase diagram. The developed model not only accounts correctly for 
the observed spin-reorientation transitions, but it predicts also their unusual angular dependence. 
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I. INTRODUCTION 

Phase diagrams of antiferromagnets in externally applied 
magnetic fields have been studied, both experimentally and 
theoretically, for more than 50 years^ In conventional 
collinear antiferromagnets the competition between Zeeman, 
anisotropy, and exchange interactions gives rise to several 
phase transitions. For example, if the antiferromagnetic or- 
der parameter has a preferred direction (easy axis), a magnetic 
field applied along that axis provokes a spin-flop transition at 
a critical field value H c , where the loss in anisotropy energy is 
compensated by a gain in Zeeman energy. Another transition 
in standard antiferromagnets is the spin-flip transition, which 
occurs when the Zeeman energy exceeds the exchange energy. 

In the more complicated cases of non-collinear and/or frus- 
trated antiferromagnets the choice of an ordered phase and 
the orientation of the relevant order parameter are dictated by 
a fine balance between the different interactions or by fluc- 
tuation effects^ This close competition implies rich phase 
diagrams with unusual features (such as, e.g., the appear- 
ance of magnetization plateaus), which have been studied both 
theoretically^ and experimentally^ 

In view of the above, it came as a surprise when the quasi- 
one-dimensional Heisenberg antiferromagnet BaCu2Si2C>7 
(hereafter BCSO), identified at low fields as an easy-axis 
collinear antiferromagnet, revealed "extra" spin-reorientation 
transitions, both in an applied field along the easy axis^ as 
well as in transverse applied fields! 10,11 By now, the phase 
diagram for applied fields along the main directions of the or- 
thorhombic crystal unit cell is well established: 12 at 2 K, with 
the field applied along the easy axis c, two spin-reorientations 
are observed (at H ci = 18.9 kOe and H c2 = 47.1 kOe), with 
the field applied along the b axis, one spin-reorientation is 
observed at H c ^ = 73.9 kOe and, finally, with the field ap- 
plied along the a direction another spin-reorientation transi- 
tion takes place at H C 4 = 114 kOe. 

On the basis of neutron scattering experiments^ the tran- 
sitions occurring in a magnetic field applied along the easy 



axis were interpreted as consecutive rotations of the order pa- 
rameter away from the easy axis, to a plane normal to it, fol- 
lowed by a rotation within this plane. A weak noncollinearity 
in the spin-flopped phase, at H c \ < H < was suggested as 
welli^ Subsequent analyses of magnetic resonance data con- 
firmed that transitions in the transverse direction were indeed 
rotations of the sublattice magnetization away from the easy 
axis. 

The temperature dependence of the static magnetization 
(see e.g. Ref.[H) exhibits unusual features at low temperatures. 
The results of neutron-scattering experiments confirmed the 
one-dimensional character of the spin system of BCSOi^ and 
the magnetic susceptibilities along the a, b and c direction 
indeed exhibit the Bonner-Fisher maxima expected for spin 
chains at elevated temperatures. 8 At lower temperatures, how- 
ever, unexpected increases of %b an d Xc with decreasing tem- 
perature are observed above the Neel temperature 7^, atyp- 
ical for this type of spin systems. More recent experiments 
probing the 29 Si NMR line shift and its temperature depen- 
dence also revealed deviations from the expected conventional 
Bonner-Fisher behaviour^ 

The theoretical description of the physics underlying these 
phase transitions is neither complete nor satisfactory. A 
phenomenological approach was used to describe the low- 
temperature phase transitions and the antiferromagnetic res- 
onance spectra. 11 ' 12 The peculiarities of the susceptibility 
above the Neel temperature were discussed in relation with 
the known behavior of weak ferromagnets4£ Yet, these ap- 
proaches predict an unusual increase of certain parameters 
with respect to their conventional estimates, hence requiring 
an additional refinement of the theory. Mean-field theory 
models have been attempted in the pastji& but they require the 
exact knowledge of many (often unavailable) microscopic pa- 
rameters. By contrast, a thermodynamics-based approach has 
better chances to capture the overall physical picture, while 
being less demanding in terms of parameter knowledge. 

Recent NMR studies — have provided a direct access to the 
local magnetization, both above and below the Neel temper- 
ature, revealing that a field-induced transverse magnetization 
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appears on the magnetic ions. The related staggered magnetic 
field is an additional parameter that needs to be considered in 
a comprehensive discussion of the low-temperature magnetic 
features ofBCSO. 

In this work, based on a conventional Ginzburg-Landau 
(GL) approach which takes into account the new experimental 
findings, we reconsider the interpretation of the phase transi- 
tions in BCSO. The field-induced transverse staggered mag- 
netization (TSM) competes with the order parameter of the 
phase with spontaneously broken symmetry below 7p>j. This 
competition may cause additional phase transitions in non- 
zero external magnetic fields and hence influence the mag- 
netization process. We demonstrate that a GL-type analysis 
of the available data provides a semi-quantitative description 
of the phase diagram and predicts both the phase boundaries 
and their variation upon changing the external magnetic-field 
orientation with respect to the crystal axes. 

In Sec. [TT] we describe and justify the phenomenological 
model, whose parameters are obtained from fitting the model 
calculations to experimental data. This is shown and dis- 
cussed in detail in Sec.|IIIJ leading finally to some conclusions 
at the end of the manuscript. 



TABLE I. Definitions of the relevant magnetic vectors and effect of 
the symmetry operations on their components. Ions positions are 
enumerated following Figure Q] Signs in the first column show 



effect of the symmetry operations / ((x,y,z) - 

({x,y,z) -> (5 -x,-y,z+z)) and C] ({x,y,z) 
correspondingly. 



{-x,j+y, -z)), 



L] = Si — S2 — S3 + S4 + S5 — S(, — S7 + Sg 
L2 = Si — S2 + S3 — S4 + S5 — S6 + S7 — Sg 

L4 = Si — S2 — S3 + S4 — S5 + Sg + S7 — Sg 

L6 = Si — S2+S3 — S4 — S5 + S6 — S7+S8 

Li*, L 2y 
L2x, Li y , H z 

Hy, L\ z 
L(,X> L4y 

L 6z 



II. THE THEORETICAL MODEL 

We start our discussion from the paramagnetic phase. As 
the temperature decreases towards the transition temperature, 
the thermodynamic functions can be expanded over powers of 
the different order parameters (i.e., over different irreducible 
representations). The symmetry group of BaCu2Si207 {Pnma 
or T>\\, see Ref. \\% includes only one-dimensional represen- 
tations which have been classified earlier i 15 i 18 Here, for the 
sake of consistency, we will use the same classification and 
axes notation (x \\ a, y \\ b, and z || c). 
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FIG. 1. Positions of the magnetic ions and their oxigen surroundings 
in the unit cell of BaCu2Si207. The fractional cell coordinates of 
the Cu 2 + ion (1) are givenil by (1/4-0.028,0.004,3/4 + 0.044). 
Chains direction is shown by bold line. 

Because of the quasi-one dimensionality of BaCu2Si207, 
its magnetic order, involving spin-i/2 Cu 2+ ions, appears 



at temperatures much lower than those implied by the in- 
chain exchange interaction strength (7n « 9.2 K"A^ while 
J = 24.1 meV— ). Thus, at the phase transition, representa- 
tions corresponding to the ferromagnetic in-chain order can 
be totally ruled out. Among the eight irreducible representa- 
tions, only four correspond to the antiferromagnetic in-chain 
ordering. These differ by the mutual orientation of the spins 
in the transverse direction and, following Ref. [l_8| we denote 
them by Li, L2, L4, and L,^. These representations and effect 
of symmetry operations on the corresponding components of 
magnetic vectors are shown in the Table U Lg represents the 
ordering pattern in the form of a ferromagnetic alignment of 
neighbouring spins along the a axis and an antiferromagnetic 
(AFM) alignment of neighbouring spins along the b axis. Li 
instead represents an AFM alignment along both the a- and 
the b axis. An AFM alignment along a, but FM alignment 
along b is represented by L2. Finally, L4 represents an FM 
alignment along both the a- and the b axis. All of them are 
consistent with an AFM order along the c axis. 

In the absence of an applied field the main order parameter 
Lg develops at the Neel point, as established by neutron scat- 
tering experiments. - Bilinear invariants that couple different 
representations with the principal order parameter and mag- 
netic field can be directly deduced from Table HJ 

LfaLfy, LfyL4 X , L\ y H z , L\ z H y , Lz x H z , L2 Z H X (1) 

Thus, while the components of L4 can appear only simulta- 
neously with the corresponding components of the main order 
parameter Lg, and are exactly zero above T^, the components 
of L] and L2 can be induced by an applied magnetic field, 
independent of the main-order parameter's existence or orien- 
tation. If Li (or L2) would be the principal order parameter, 
invariants L\ y H z and L\ z H y (or L2 X H Z and L2 Z H X ) would lead 
to the appearance of the weak ferromagnetism, as it happens 
in the related compound BaCu2Ge207^ 
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The representations Li, h 2 , and L4 do not need to be in- 
cluded in the free-energy expansion if only the macroscopic 
energy is of interest, minimization over their components re- 
sults in the thermodynamic function expansion over principal 
order parameter. At the microscopic level, however, these rep- 
resentations are components of the local magnetization which 
are experimentally accessible. 

In particular, recent NMR studies^ have demonstrated the 
presence of a field-induced staggered magnetization, both 
above and below 7n. A similar effect is well known in the case 
of weak ferromagnets above the transition temperature! 20 ' 21 In 
order to capture this situation in our model, the field-induced 
transverse staggered magnetization, here represented by Li 
and L2, have to be included in the GL free-energy expansion. 
For our purposes and in order to simplify the calculations, we 
refrain from considering L4, which would only lead to a slight 
renormalization of the anisotropy constants. 

The value of the critical exponent describing the growth of 
L6 close to Tn, as determined by neutron scattering-^ and by 
NMRii is close to 0.25, differing significantly from the clas- 
sical GL value of 0.5. Therefore, a GL approach in dealing 
with the present case has its limitations. Yet, it can still pro- 
vide useful insights into the nature of the phase transitions and 
the understanding of competing magnetic order parameters of 
quasi-one-dimensional spin systems. 

Following the general theory^ we use the thermodynamic 
function <t>, defined in such a way that <9<t>/<9H = — H/(4n) — 
M: 
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The quadratic terms A,L? describe the exchange rigidity to- 
wards the formation of the corresponding order parameter. As 

usual, Ag = ^(T — 7^°'), while A\ and A 2 remain positive. 
We suppose that, because of the low dimensionality of the spin 
system, A 1.2 are particulary small close to Their temper- 
ature dependence can be relatively strong, however, and thus 
we assume a linear temperature dependence in the vicinity of 

T N , i.e.: A L2 (T) = A§ [1 + Zff(T - T< 0) )]. All other co- 
efficients are postulated to be temperature independent. The 
fourth-order term B^\J^ fixes the magnitude of the principal 
order parameter Lg below the transition. The following terms 
B,gLgL? and B' j6 (L^ -L,) 2 describe the competition between 
the field-driven TSM Li 2 and the exchange-driven principal 
order parameter Lg, the key topic of our paper. The next 
terms describe the usual exchange contributions to the mag- 
netization (D and D'), the anisotropic interactions affecting 
the principal order parameters (a x and a y ), and the anisotropic 



interactions responsible for the coupling of Li ,2 to the mag- 
netic field (a,, z and p x ,z)- Finally, the last line represents the 
paramagnetic susceptibility % p and its corrections (e.g., due to 
g-factor anisotropy) y x and J y . 

The magnetic anisotropy axes of BaCu2Si207 have been 
identified by magnetization,— neutron scattering^ and antifer- 
romagnetic resonance (AFMR) 11 ' 22 experiments: the c axis 
represents the easy axis, while b is the secondary easy axis 
(hence implying a x > a y > 0). The B,g coefficients are ex- 
pected to be positive when non-coexisting representations are 
favored. The B' j6 coefficients instead, which are responsible 
for the preferred mutual orientation of the principal order pa- 
rameter and TSM, do not have an a priori given sign. 

The field-induced TSMs L12 can be found by minimizing 
the value of 4>, as given by Eq. ©. For instance, when H || y 
the only non-zero component is 



a y H y 



2(A l+ B l6 Ll+B\ 6 Ll) 



(3) 



This TSM can be induced by the magnetic field already above 
Tn. At the Neel point Lg starts to grow and to suppress L\ z 
( A 1 , Z? 1 6 > 0). Similarly, for the other two principal field ori- 
entations, H || x induces L 2z , while H || z induces both L\ y and 
L 2x . All the induced transverse staggered magnetizations de- 
pend linearly on the applied field. 

The substitution of the field-induced TSMs found above 
into Eq. (O results in field-dependent terms quadratic in H. 
These terms correspond to the corrections to the susceptibility 
which, in case of applied fields along the crystalline axes, can 
be written as: 
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(4) 
(5) 



(6) 



These corrections to the susceptibility provide a natural ex- 
planation for the additional contributions to % observed above 
Tjsr. Below Tn the Ax terms depend on the orientation of the 
main order parameter and, hence, provide clues about the pos- 
sible spin-reorientation transitions. For instance, for H || x, a 
transition with the rotation of the order parameter from the 
easy-axis z towards the v-axis is possible only if B' 26 > 0. The 
corresponding transition field is 



H, 



\a y A\ 



c4 



(7) 



2h 



Note that the a y and j3 2 parameters are of the same order of 
magnitude as the corrections due to the spin-orbit interaction. 
Consequently, in the ordinary antiferromagnet, the field H C 4 
should be comparable with the exchange field. The relatively 
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small (as compared with the exchange field defined by the in- 
chain exchange integral) value of H C 4 is in fact due to the tiny 
value of A2, in turn related to the one-dimensionality of the 
system. The positiveness of B' 26 means that at high applied 
fields an orthogonal alignment of the field-induced TSM L2 
and of the exchange-driven Lg order parameters is favored, 
while at low fields (for H || x) both L2 and Lg are parallel 
to z. As a result, since the field-induced order parameter is 
determined by the applied field, the main order parameter will 
start to rotate whenever the susceptibility -related energy gains 
overcome the anisotropy-related losses. 

Similarly, for H || y, B' l6 > will cause a spin-reorientation 
at a critical field 



H, 



c3 



Aa x A\ 



(8) 
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exchange term D'H 2 , which again can be explained by the 
particular smallness of the A\ 2 parameters. 

The paramagnetic-antiferromagnetic phase boundary was 
studied in the mean-field approach considering coupled chains 
in a staggered fieldJ£ This approach demonstrated the sup- 
pression of the symmetry-breaking ordered phase by the stag- 
gered field. Our results show a similar behavior: positive 
Bis and B26 (which corresponds to the competition between 
the principal order parameter and the TSM) leads to the de- 
crease of the Neel temperature with respect to the ordinary 
case Z?i6 26 = 0. Besides that, our results show that under spe- 
cific conditions (corresponding to the positiveness of B\ 6 26 m 
our thermodynamic model) the transverse staggered field af- 
fects not only the stability of the ordered phase but necessar- 
ily leads to new spin-reorientation transitions in the ordered 
phase. 



Finally, for H || z, a rotation of the main order parameter from 
the secondary easy-axis y towards the hard axis x is possible 
if 
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Cl ^B'^Al-^AY 
The normal spin-flop occurs at the field: 
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Likewise, field-induced shifts of the Neel temperature can 
be straightforwardly calculated from Eq. (0. As an example, 
we consider the main order parameter to be oriented as in the 
case of BaCu2Si2C>7 and obtain: 
For H || x and H > H c4 (i.e., L 6 || y): 



III. COMPARISON OF THE MODEL WITH 
EXPERIMENTAL DATA AND DISCUSSION 

^ To compare the theoretical model with experimental data 
we will rely mostly on the already published phase diagram 12 
and on NMR and magnetization data^i 4 - 

The expansion (O of the thermodynamic function <t> in- 

(10) eludes 19 explicit parameters and 2 coefficients ^ r 2 \ which 
describe the temperature dependence of A \ 2 . Although the 
existing experimental data allow us to fix all the parameters, 
part of them, however, are not critical for the computation of 
the phase diagram. Since our procedure is affected by the 
choice of certain extrapolations (see below), the values used 
here differ slightly from those reported in Ref . (3- 

Since all the measurable parameters, except the absolute 
magnitudes of the field-induced order parameters, depend 
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For H || y and H > H c3 (i.e., L 6 || x): 
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For H || z and H > H c2 (i.e., L 6 || x): 
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In case of BaCu2Si2C>7 the Neel temperature was found to 
increase when H || x and to decrease for applied fields along 
the other two directions. This means that additional correc- 
tions turn out to be comparable in magnitude with the main 



only on the ratios a 2 z /Ai, j3 2 ./A 2 , B16/A1, B\JA\, B 2 (,/A 2 , 

and B' 26 /A 2 , the values of the A\ 2 terms at 7n (a\ l) can be 
evaluated independently. 

By requiring that the saturated value of the main order pa- 
rameter at zero applied field is unity and by using the known 
magnitude of the zero-field specific heat jump^ AC = 0.61 
J/(K-mol) (per mole of compound), we obtain = 0.305 x 
10 7 emu/(K- mol Cu) and B 6 = 1.40 x 10 7 emu/(mol Cu). 
Paramagnetic contributions to the susceptibility can be de- 
termined by an extrapolation of the susceptibility curve for 
ID Heisenberg chains.— To this end we used the exchange 
integral value J = 24.1 meV, determined by neutron scatter- 
ing experiments, 13 and fitted the high-temperature tails of the 
measured %(T) curves by setting the g-factor values equal 
to 2.20, 2.00 and 2.06 for H || a, b and c, respectively (see 
Fig. |2]i. This extrapolation at the Neel temperature yields 
X P = 6.29 x 10~ 4 emu/(mol Cu), y x = 0.88 x 10~ 4 emu/(mol 
Cu), Yy = -0.34 x 10~ 4 emu/(mol Cu). 

A subtraction of the corresponding Bonner-Fisher curve 
from the measured magnetization data (see Fig. |2]i provides 
the additional contributions to the susceptibility, as described 
by Eqs. ©-(O. By extrapolating the latter curves to the 
Tj<j value, we can fix the following parameters combinations: 
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FIG. 2. (color online) Upper panel: Temperature dependence of the 
magnetic susceptibilities measured at 10 kOe (Ref. [l4h (symbols). 
Model curves for a pure ID Heisenberg antiferromagnet (Ref. l23h 
with J = 24.1 meV, g a = 2.20, g b = 2.00, and g c = 2.06 (dashed 
lines). Lower panel: Difference between the measured %(T) and the 
corresponding model curves (open symbols), variation of the NMR 
shifti 4 . with temperature (below the average position of the split 
peaks is shown) (closed symbols), model calculations of the present 
work (solid lines). 



4fc(r N ) = &= 0.89 x 10- 4 emu/(mol Cu), A^(T N ) 
■ emu/(mol Cu) and A&(r N ) = + 



^ =2.28 x lfr 4 — 1 rwv™, i a., tizA - «! , ft 2 

5.20 x 10~ 4 emu/(mol Cu). The coefficients <§J '2 can be es- 
timated from the temperature dependence of the NMR shifts 
across the transition and from the additional contributions to 



the susceptibility above c§j = 0.037 K _1 and ^ 
0.02 K . The ratio of the anisotropy constants a x /a y = 3.4 
is known from AFMR data. 1 1 Finally, the main correction to 
the susceptibility, due to the onset of an antiferromagnetic or- 
der, D, can be estimated from the magnetization curves as 
D = 8.0 x 10~ 4 emu/(mol Cu). 

The few remaining parameters can be tuned to the extent 
that the calculated critical-field values and the field-induced 
shifts of 7n agree best with those obtained from the experi- 
ment. The critical fields at the paramagnet-antiferromagnet 
phase boundary are H c \ = 18 kOe, H C 2 = 53 kOe, H c t, = 
81.5 kOe, and H c \ =116 kOe, while the values of the Neel 
temperature shifts at 14 T (AT N = r^ 14T) - 7^ 0) are 0.16(2) K 
for H || x, -0.1 1 (4) K for H || y, and -0.98(2) K for H |j z)H 
All the parameter values and their relevant combinations are 
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TABLE II. Parameters of the thermodynamic function <I>, as given by 
Eq. J2f, used in the model calculations. 
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6.29 x 10~ 4 emu/(mol Cu) 
0.88 x 10~ 4 emu/(mol Cu) 
-0.34 x 10~ 4 emu/(mol Cu) 

8.0 x 10~ 4 emu/(mol Cu) 
-0.88 x 10~ 4 emu/(mol Cu) 
0.037 K _1 
0.02 K" 1 
0.305 x 10 7 emu/(K- mol Cu) 
1.40 x 10 7 emu/(mol Cu) 
1.78 x 10~ 4 emu/(mol Cu) 
4.56 x 10~ 4 emu/(mol Cu) 
2.03 x 10~ 4 emu/(mol Cu) 
8.37 x 10~ 4 emu/(mol Cu) 
0.622 
1.197 
0.892 
0.332 

6.76 x 10 5 emu/(mol Cu) 
1.99 x 10 5 emu/(molCu) 
1.06 x 10 7 emu/(mol Cu) 
6.63 x 10 5 emu/(mol Cu) 
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FIG. 3. (color online) Experimental phase diagram of BaCu2Si207 
(Ref. [T2I) (magnetization data — open circles and squares, specific 
heat measurements — closed triangles) vs. modelled phase bound- 
aries (solid lines). Orientations of the principal order parameter L6 
and of the field induced TSM Li 2 are shown for the corresponding 
phases. 



listed in Table|TTl 

With these parameters the magnetization curves can be cal- 
culated and the relevant phase boundaries can be established. 
Figures [2] and [3] show that the modelled curves are reason- 
ably close to experiment. The main failure of the model con- 
sists in the predicted temperature dependence of the high-field 
spin-reorientation transitions, which is not observed in the ex- 
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FIG. 4. (color online) Upper panel: Comparison of the modelled 
temperature dependence of the field-induced order parameter (solid 
line) vs. measured NMR line shifti 4 . (HqH = 7.02 T, H || b) (circles). 
Lower panel: Temperature dependence of the induced order param- 
eters for H || c and [l H = 4.0 T. 



periment. Probably this can be accounted for by consider- 
ing the role of thermal fluctuations, which are known to sta- 
bilize collinear magnetic structures 4 - (field-induced TSM and 
exchange-driven order parameter are collinear in the low-field 
phases for H || a,b). Our model provides an adequate descrip- 
tion also for the observed NMR shifts, both above and be- 
low the Neel point, which, for H || b,— is mostly due to the 
staggered magnetization pattern related to L\ z (see Fig. |4). 
The scaling of the NMR shift 14 and of A% y above 7n (see 
Fig.|2]i confirms once more the reliability of the model [com- 
pare Eqs. (O and @]. 

The evaluation of the individual values of Aj ^ mentioned 
above is possible if the values of the field-induced TSM's at 
the phase transitions in differently oriented magnetic fields are 
known. The evaluation of the latter is not straightforward but 
still feasible. 

An estimate for the field-induced TSM, L\ z = 0.035 fis per 
Cu 2+ ion (in an applied field of 7.02 T), was obtained in re- 
cent NMR work J 4 Since in our model the saturated value 
of the main order parameter is normalized, while its mea- 
sured valued is ca. 0.15 JJ-b per Cu 2+ ion, this corresponds 
to Li z (7n) = 0.23 in normalized units. From the latter value 
and from Eq. (0 we find A ( j 0) = 1.06 x 10 7 emu/(mol Cu). 

The value can be estimated from the canting of 
the magnetic structure, as observed in neutron scattering 
experiments^ with H || c. These experiments revealed an L2 X 
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FIG. 5. (color online) Modelled angular dependence of the spin- 
reorientation critical fields, performed at three different tempera- 
tures. The missing data for T = 9 K in the rotation from b- toward 
the c-axis reflect the crossing of the PM-AFM boundary before a 
spin-reorientation could occur. 



component of ca. 0.17 (corresponding to a canting angle of 
~ 10° at )1qH ~ 4 T). However, because of the coincidence of 
the structural and magnetic Bragg peaks in BaCu2Si207, the 
determination of the magnetic scattering intensities required 
the subtraction of the peak intensities recorded just above 7^. 
Since the contribution of the field-induced order parameter is 
subtracted as well, the observed low-temperature magnitude 
of L2 X represents, in fact, only the change of L^x across due 
to the competition between the field-induced and exchange- 
driven orders. Our calculations (Fig. HJi show that Li x at 6 K 
amounts to ~ 50% of its value at 7n- Therefore, by assum- 
ing L2jc(7n) = 0.35, we obtain A 2 °^ = 6.63 x 10 s emu/(mol 
Cu). Note that an L\ y component should also exist for H || c. 
However, its magnitude at 7n and HqH = 4 T is, as calcu- 
lated using the found parameters, ca. 0.17 (i.e. only half of 
L2x). Besides, our model calculations show that L\ y strongly 
changes below only at H c \ < H < H C 2 when the princi- 
pal order parameter is also aligned along the y-axis. These 
reasons probably explain why L\ y has not been observed ex- 
perimentally. The above estimates for A J ] are consistent with 
their expected small values. In fact, both of them are compa- 
rable with £ 6 (T-T N ) evaluated at (T — r N ) ~1K. 

Finally, we also computed the angular dependence of the 
critical fields for spin-reorientation transitions (see Fig.[5}. All 
the marked fields correspond to real transitions, accompanied 
by a jump in the susceptibility and by a sudden reorientation 
of the main order parameter To ensure that the angular 
dependence is not affected by changes in magnitude of the 
main order parameter Lg (due to closeness to the PM-AFM 
phase boundary), the modelling was carried out at different 
temperatures (6, 8, and 9 K). All curves are qualitatively sim- 
ilar, with differences mostly due to the above mentioned tem- 
perature dependence of the upper critical fields H c % 3 4 and, 
partially (for the 9 K curve), to the crossing of the PM-AFM 
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phase boundary before a spin-reorientation transition has oc- 
curred. Our model predicts quite a remarkable angular depen- 
dence for the critical fields. On rotating from the hard axis a 
towards the secondary easy-axis b, the critical field H c 4 trans- 
forms smoothly into H c i . Upon further rotation towards the 
easy axis c, the field required for a spin-reorientation transi- 
tion first grows very rapidly and, at an intermediate critical 
angle, diverges to infinity. Subsequently, it reappears from 
the high-field zone and finally converges to the critical field 
H C 2- A rotation from the easy axis c towards the hard axis a 
(right panel in Fig. |5]l demonstrates the merging of the two 
critical fields H c \ and H C 2 at a certain angle. Then the spin- 
reorientation transition reappears from high fields and, finally, 
close to the a axis, evolves towards H c 4. To the best of our 
knowledge the angular dependence of the critical fields in 
BaCu2Si207 has not yet been studied. Its experimental in- 
vestigation would represent an independent additional test of 
the proposed model. 



IV. CONCLUSIONS 

By making use of the Ginzburg-Landau theory of phase 
transitions we propose a semi-quantitative description of the 
magnetic phase diagram of ID systems with competing inter- 
actions. We have demonstrated that the competition between 
the field-induced and the exchange-driven order parameters in 
a quasi-one-dimensional antiferromagnet can lead to an un- 
usual phase diagram and to remarkable deviations of the mag- 
netization process from that expected in a ID Heisenberg an- 
tiferromagnet. Additionally, in the BaCu2Si207 model sys- 
tem, we predict an unusual angular dependence of the critical 
fields, which will be object of future experimental investiga- 
tions. 
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